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The magnetic excitations of the iron pnictides are explained within a degenerate double-exchange 
model. The local-moment spins are coupled by superexchanges Ji and J2 between nearest and 
next-nearest neighbors, respectively, and interact with the itinerant electrons of the degenerate d^z 
and dyz orbitals via a ferromagnetic Hund exchange. The latter stabilizes (7r,Q) stripe antiferro- 
magnetism due to emergent ferro-orbital order and the resulting kinetic energy gain by hopping 
preferably along the ferromagnetic spin direction. Taking the quantum nature of the spins into 
account, we calculate the magnetic excitation spectra in the presence of both, super- and double- 
exchange. A dramatic increase of the spin-wave energies at the competing Neel ordering wave vector 
is found, in agreement with recent neutron scattering data. The spectra are fitted to a spin-only 
model with a strong spatial anisotropy and additional longer ranged couplings along the ferromag- 
netic chains. Over a realistic parameter range, the effective couplings along the chains are negative 
corresponding to unfrustrated stripe antiferromagnetism. 
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I. INTRODUCTION 

The discovery of superconductivity in the pnictides^^^ 
with transition temperatures challenging those of single- 
layer, high-Tc cuprates, immediately raised the ques- 
tion of whether, despite all their differences, the two 
classes of materials share the same key mechanism for 
superconductivityji^— Arguably the most striking simi- 
larity is that superconductivity emerges upon doping an- 
tiferromagnetically ordered parent compounds. In the 
pnictides, however, the magnetic ordering is unusual, 
characterized by an antiferromagnetic arrangement of 
ferromagnetic chains, corresponding to an in-plane or- 
dering wave vector Q = (tt, 0).Jfi"i8. Whereas the pnic- 
tides are metallic, the parent cuprates are Mott insula- 
tors. The conductivities in the pnictides are typical of 
bad metals, suggesting that electronic correlations^^TiSS 
are crucial. By contrast, the local-density approximation 
seems to be adequate^S to describe their band structure. 

As a result of this dichotomy, both itinerant-magne- 
tisni^Sr^ and local-moment^"— i^Sti^ scenarios have 
been suggested to explain the unusual stripe antiferro- 
magnetism. Although the former weak-coupling theories 
which attribute the magnetism to a spin-density- wave in- 
stability of a nested Fermi surface, can explain the mag- 
netic low-energy excitations around the correct ordering 
wave vector, they fail to describe the spectra at higher 
energies which have been measured in great detail by 
inelastic neutron scattering 4^>^ In particular, the itin- 
erant scenarios can not explain the observed maximum of 
the dispersion at (tt, tt) but rather suggest that the exci- 
tations rapidly dissolve into a particle-hole continuum22 
which has not been found experimentally up to energies 
of 200meV. So far, a consistent description of the excita- 
tions over the entire energy range has been obtained only 
by using suitably parametrized Heisenberg models. 

Because of the positions of the arsenic ions above 



or below the iron plaquettes, such a spin-only model 
is expected to be strongly frustrated with compara- 
ble, antiferromagnetic superexchanges Ji and J2 be- 
tween nearest and next-nearest neighbors. In this 
regime, the model indeed exhibits long-range stripe- 
antiferromagnetic order^"*^"*^ and the strong frustration 
and proximity to a continuous magnetic phase transition 
might explain why the observed magnetic moments are 
relatively small. Interestingly, the neutron-scattering 
experiments tell a radically different story. The spin- 
wave velocities indicate that the exchange coupling along 
the ferromagnetic spin direction is much smaller than the 
one perpendicular to the chains?^ Jiy <^ Jix- More re- 
cently, it has been argued^ that the observed maximum 
of the dispersion at (tt, tt) requires an even slightly ferro- 
magnetic exchange Jiy < corresponding to an unfrus- 
trated spin model in contrast to early claims^^ of high 
frustration. 

What might be the cause of such a strong spatial 
anisotropy? In fact, before the magnetic order sets in, 
a structural transitions occurs at which the two in-plane 
lattice constants become inequivalent. The structural 
and magnetic transition are clearly separated in the so- 
called '1111 compounds' fiS^J^ whereas they occur at the 
same temperature in the '122 family'.— "= However, in- 
specting the numbers, it appears that the orthorhombic 
lattice distortion is too small, by two orders of magni- 
tude, to explain the magnetic anisotropy<^ 

To this end, some have proposed^i^"— that orbital- 
ordering physics of a similar kind as in the manganite 
transition-metal oxides not only provides a mechanism 
for the lattice distortion but more importantly explains 
the strong in-plane anisotropics. In particular, it has 
been argued^ that due to an orbital degeneracy, the lo- 
calized limit is described by a complicated spin-orbital 
superexchange (Kugel-Khomskii) model rather than by a 
Heisenberg Hamiltonian. Further, it was shown that the 
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stripe antiferromagnetism is stabilized by ferro-orbital or- 
der which breaks the in-plane lattice symmetry and in- 
duces a strong anisotropy between the magnetic exchange 
couplings. 

Since the C4 lattice symmetry is broken by the or- 
bital order and the accompanying orthorhombic dis- 
tortion, the electronic structure must reflect this spa- 
tial anisotropy with reduced symmetry. Indeed, such 
an anisotropic electronic state has been confirmed re- 
cently in scanning-tunneling-microscopy (STM)^* and 
in-plane resistivity measurements . 1 Dramatic Fermi- 
surface reconstructions at the structural transition^ as 
well as enormous transport^ and phonon^ anomalies, 
have been interpreted as indirect evidence for orbital or- 
dering. Hence, on this interpretation, it is the orbital 
ordering that underlies the electronic anisotropy, not an 
inherent anisotropy due entirely to the electrons, indica- 
tive of a true nematic state. 

Consequently, an open problem with the pnictides 
is the role itineracy and local physics play in mediat- 
ing the apparently unfrustrated anisotropic magnetism. 
In this work, we start from the idea of the 'local- 
itinerant dichotomy —i^i^ of the iron pnictides and mo- 
tivate an effective degenerate double-exchange (DDEX) 
model, very similar to the ones used to describe metallic 
manganites with orbital degeneracies 1^"— To be more 
precise, we assume a ferromagnetic Hund coupling be- 
tween the itinerant bands of the doubly-degenerate d^z 
and dyz orbitals and the local moments, which are de- 
scribed by the aforementioned J1-J2 Heisenberg model. 
In the context of the manganites it has been shown^"— 
that such DDEX models exhibit phases with long-range 
stripe-antiferromagnetic order. Despite the antiferro- 
magnetic couplings between local-moment spins, ferro- 
magnetic spin chains are stabilized by emergent ferro- 
orbital order. In this phase, the itinerant electrons are 
directed predominantly along the chains which minimizes 
the kinetic energy and gives rise to a highly anisotropic 
electronic state. Moreover, the double exchange is ex- 
pected to strongly suppress the effective coupling be- 
tween local moments along the chains, and possibly to 
make it even ferromagnetic. 

In this work, we analyze the effective DDEX model 
and indeed find that the orbitally ordered (tt, 0) antifer- 
romagnet is stable over a realistic parameter range for the 
parent iron-pnictide materials. In particular, the seizable 
next-nearest neighbor superexchange J2 further stabilizes 
this phase. Whereas these results are to a large extent 
not surprising given the similarities with the manganites, 
the magnetic excitation spectra so far have been calcu- 
lated only for a ladder system;^ For the manganites, the 
DDEX model is usually simplifiedSSr— by treating the 
core spins as classical and by taking the limit of large 
or infinite Hund's coupling Jh which is not justified for 
the pnictides. Here, we instead focus on the regime of 
small and intermediate Jh and develop the tools to cal- 
culate the magnetic excitation spectra in the presence 
of both, super- and double-exchange to linear-spin-wave 



order. The spectra are found to be in good agreement 
with the neutron scattering data4^ In particular, in some 
parameter space the double exchange along the ferromag- 
netic chains can overcompensate the bare antiferromag- 
netic superexchange as suggested by the experiment. 

This paper is organized as follows. In Sec. [Hi we con- 
struct the local- itinerant, DDEX model. Sec. IIIII deals 
with the methods we use to calculate the magnetic exci- 
tation spectra. In Sec. IIVI we summarize our results, in- 
cluding ferro-orbital ordering, the spin-wave dispersions, 
and the magnetic anisotropics. Finally Sec. |V] discusses 
several aspects of our theory and validates it's applica- 
bility to the pnictides. 



II. MODEL 

In this Section, we proceed to motivate the DDEX 
model for the pnictides. This model accounts for the 
presence of local moments, as suggested by the neutron- 
scattering experiments, as well as itinerant electrons re- 
sponsible for the bad-metal behavior of the parent com- 
pounds. Moreover, the orbital degeneracy in combina- 
tion with Hund's coupling between electronic and spin 
degrees of freedom gives rise to orbital-ordering physics 
beyond simple band-structure theory. The Hamiltonian 
consists of three parts, 

n = nioc + nit + nH, (i) 

where Hioc describes the superexchange couplings be- 
tween local moments, T-Lu the itinerant electrons of the 
degenerate dxz and dyz orbitals, and Hh the ferromag- 
netic Hund coupling between local moments and itiner- 
ant electrons. In order for this model to be valid, Hund's 
coupling Jh should be small compared to the tetrahedral 
crystal-field splitting between the t2g and Cg multiplets, 
but larger than the tetragonal splitting between the dxy 
orbital and the degenerate d^zAyz doublet.^- 

The local moments with spin S are coupled by su- 
perexchanges Ji and J2 between nearest and next-nearest 
neighbors, respectively. The corresponding Heisenberg 
Hamiltonian reads 

where, for convenience, the superexchanges are measured 
in units of S^. Likewise, Hund's exchange Jh which cou- 
ples the electron spins to the local moments will be mea- 
sured in units of S. This convention will facilitate our 
large- S* expansion later. The superexchanges are medi- 
ated by virtual hopping processes via the p-orbitals of the 
arsenic ions which have alternating positions above and 
below the iron plaquettes. Certainly, a quantitative com- 
parison of the exchange couplings would require knowl- 
edge of the two different Fc-As-Fe bond angles and the 
precise shape of the orbitals. Assuming that the virtual 
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FIG. 1; (Color online) Illustration of the degenerate double- 
exchange model for the pnictides. (a) The local moments are 
coupled by nearest- and next-nearest-neighbor exchanges Ji 
and J2, respectively, and interact via a ferromagnetic Hund 
coupling Jh with the itinerant electrons of the degenerate 
dxz, dyz orbitals. (b) Illustration of the hopping parameters 
in a two-band model of these orbitals (shown as projections in 
the plane), (c) Resulting ferro-orbital order which stabilizes 
(tt, 0) antiferromagnetism by directing the kinetic energy of 
the electrons along the ferromagnetic spin direction. 



processes for the nearest-neighbor and diagonal bonds in- 
volve roughly the same energies, we estimate Ji ~ 2J2 
since two exchange paths via different arsenic ions con- 
tribute to Ji . Therefore, we expect the Heisenberg model 
([2]) to be strongly frustrated and potentially in the quan- 
tum disordered regime. 

The itinerant electrons of the degenerate dxz and dyz 
orbitals are described by a tight-binding Hamiltonian 

"Wit = - ^"/4aL.Cj/3,„ (3) 

where cj^^ creates an electron with spin v at site i on 
orbital a. The hopping integrals t"^^ are illustrated in 
Fig. W[h) and defined in the same way as in Ref. IH. For 
simplicity, we neglect inter- and intra-orbital Coulomb 
repulsions^ between the itinerant electrons. This is jus- 
tified since on the level of the effective DDEX model, 
the local moments are formed as a consequence of strong 
correlations whereas the residual charge carriers should 



be viewed as weakly interacting quasiparticles. We will 
assume ti to be the dominant hopping because a larger 
wave- function overlap is expected when the orbitals point 
towards one another. However, the precise shape of the 
orbitals is not determined by geometry but depends on 
quantum chemistry, for example, on the amount of hy- 
bridization between the Fe d- and As p-orbitals.^° Here 
we simply denote the orbitals by dxz and dyz because of 
the spatial symmetry shared with the atomic Fe orbitals. 
Recently, it has been suggested^ that the hybridization 
leads to a strong deformation of the orbitals which make 
the TT overlap, t2 the dominant one. We point out that 
our results persist for the exchange of ti and t2, the only 
difference being that the orbital polarization will be in- 
verted in order to maximize the overlap along the ferro- 
magnetic spin direction. 

We do not attempt to fit our hopping parameters to 
reproduce the electron and hole pockets as has been done 
in previous two-orbital models^i^ at the level of a tight- 
binding approximation. Such parameters are inherently 
arbitrary since the hopping amplitudes are not uniquely 
determined by a particular constant energy culi^ and 
because Coulomb interactions and Hund's exchange are 
of the order of the electronic bandwidthi^"— Moreover, 
the Fermi surfaces in the antiferromagnetically ordered 
phase have not been clearly established yet. 

Finally the local moments and the itinerant electrons 
interact by a ferromagnetic Hund coupling, 

"^H = S^■ cl^^Cr^^'C^a^i, (4) 

where cr^y/ = (cr^, cr^, cr^)iyi/' with cr" the standard Pauli 
matrices and Jh > 0. As mentioned before, Hund's ex- 
change is measured in units of the local moment S. We 
note that similar models have been proposed in the con- 
text of the pnictides i^i^ However, the orbital degener- 
acy which is the pre-requisite for orbital-ordering physics 
and the resulting spatial anisotropics has not been in- 
cluded. 



III. METHOD 

In this Section, we outline the approximations and 
transformation we employ to analyze the complicated 
DDEX model for the pnictides. In similar models for the 
maganites, the problem is typically simplified by treat- 
ing the local moments as classical spins and assuming an 
infinitely strong Hund exchange In the pnictides, 
these approximations are not justified since the Hund 
coupling is of the order of the electronic bandwidth and 
since the local moments are small and presumably best 
described by the extreme quantum limit, S ~ 1/2. More- 
over, our goal is the calculation of the magnetic excita- 
tion spectra which will require the inclusion of quantum 
fiuctuations of the spins. Although the spins are small 
and the Heisenberg model "Hioc is strongly frustrated, it 
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is reasonable to treat the local moments on the level of 
linear spin- wave theory since the double exchange is ex- 
pected to lead to a dramatic stabilization of the magnetic 
order. Moreover, it has been argued that the 1/S expan- 
sion is much better behaved for (tt, 0) order as compared 
to conventional (tt, tt) Neel antiferromagnetism4^ 

Since Hund's coupling "Hh does not conserve the 
magnons describing the spin-wave excitation of the iso- 
lated local moments, we perform a canonical transforma- 
tion in order to identify the true magnons of the coupled 
system. Readers not interested in the details of this cal- 
culation can skip immediately to the results. Sec. IIVI 



A. Operator rotations 

Following the standard treatment of the antifcrromag- 
netic Heisenberg model, we perform the spin rotation 
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KiS-, and Sf 



KiSf where 



exp {iQ ■ Ti) = ±1 for sublattices A and B, respectively 
[see Fig. [TJa)]. Representing the rotated spin operators 
Si by Holstein-Primakoff (HP) bosons a^, a|, to the lead- 



ing order, = S- 



- a J a,- 



2Sai, and S, 



2S 



(Sf = SfztiSf), we immediately derive the following ex- 
pression for "Hioc in the linear spin-wave approximation, 

^io"c = E [^o(9)(4«? + «-9a-g) 
q 

+Baiq) (alalg + a^gag'j , (5) 

where Ao{q) = {Jicosqy + 2J2)/5 and Bo{q) = 
( Ji cos qx + 2J2 cos qx cos qy)/S. 

In order to leave Hund's coupling term Hjj invariant, 
we perform exactly the same rotation of the electron spins 
Sta = jY^uv' '^laij'^v'^'^^av' ■ This is achieved by trans- 
forming the fermion operators as ciav = Ciau for sites i 
on sublattice A and Ciav = Ciau on sublattice B. In the 
latter expression, we have defined v for u and vice 
versa. In terms of the HP boson creation and annihila- 
tion operators a^, a\ and rotated fermion operators Ciav^ 



Hund's exchange can be written as 
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where u = ±1 for up and down spins, respectively. Note 
that 'Hy ^ only involves the electronic operators and rep- 
resents the zeroth-order corrections to the electron en- 
ergies by the classical background stripe antiferromag- 



netism. "H^^ and H)^' are the couplings between the 
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electrons and the HP bosons, linear and quadratic in the 
boson operators. In the following, we include the term 
h'^^ in the itinerant-electron Hamiltonian, yielding the 
effective free-electron Hamiltonian 



He — 



17 ^ 'X/^^^ 

tit + rtjj 
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yz for a = xz and vice versa. We have defined 
—2t2COsky, s\^{k) = —2tiC0sky, £|^(fc) — 



where a = 

£f^(fc) = 

—2ticoskx — 4^3 cos fca; cos fcj, , e|^(fc) — —2^2 cos fc^; 
4^3 COS kx COS ky, and Ssik) — — ^4 sin kx sin ky. 



B. Canonical transformation 

Apparently, the interaction term h'^^ is linear in the 
HP-boson operators, which shows that these bosons do 
not represent the Goldstone modes of the system, namely 
the transverse fiuctuations of the total staggered mag- 
netic moments, which consist of not only the local mo- 
ments, but also the spins of the itinerant electrons. In 
order to correctly identify the true magnons and carry 
out the subsequent spin-wave calculations, we need to 
perform a canonical transformation of the original Hamil- 



tonian H 



-T/SW 

''-loc 



He + H. 



n 



(2) 
H ' 



n + [A,n] + -[A,[A,n]] 



(11) 



with A a suitable anti-Hermitian operator, = —A, 
such that in the transformed "H', the terms linear in a^'s 
are eliminated. Similar canonical transformations^^iZS 
and equivalent perturbative methods^i have been used to 
explain ferromagnetism in double-exchange models with 
a single itinerant band. Up to the leading order, the 
transformation is determined by 



[A, He 



(12) 



To find the A satisfying , we first diagonalize He by 
a unitary transformation Ckav = J2n Uavik)dnk, yielding 
He — k^n{k)d\^dnk- Hcrc, n labels the four elec- 
tronic bands arising after diagonalization from the two 
orbital and two spin degrees of freedom. In the new ba- 
sis of dnk^ it is easy to verify that (IT^ is solved by 



'^q'^m.k+q'^rik 



25,^ \En{k)~E^{k + q) 



X u^xik- 



g)C/^(fc) - h.c. 



(13) 
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After the canonical transformation (|TT|) . the Hamilto- 
nian up to order l/S reads U' =U^+ Hf^^ + u'i^ + U'^, 
where n'^ = [^Ml^] + i[A[A,Hc]] = i[A,-H[j^^]. The 

commutators [A,?^;*^] and [AjH^"^] are of higher orders 
in 1/5* and the boson operators, and thus can be dropped 
in the hnear spin-wave approximation. The contribu- 
tions and 'H'2 are bilinear in both, the bosonic and 
fermionic operators. After taking the expectation values 
of the electronic operators with respect to the diagonal 
electronic Hamiltonian "He, we obtain the final spin- wave 
Hamiltonian 
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(2) 
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(1)- 



E 

9 



alttq 



-B{q) [ala\ 



(14) 



with A{q) = Ao{q) + Ai + A2{q) and B{q) = 
Bo{q) + B2{q)- The constant 'self-energy' correction, 
Ai arises from whereas = i[A,'Hy^] gener- 

ates momentum-dependent corrections to both, the 'self- 
energy' and the 'anomalous amplitude', which are de- 
noted as A2{q) and i?2(<z)- These corrections are ex- 
pressed as 



Ai = 



A2{q) 



k.n 

Jn fnjk) - fmjk + q) 
2S E„ik) - Emik + q) 



(15) 



(16) 



B2{q) 



fn{k)-f„r{k + q) 



(17) 
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\ " 

2S ^ Er,{k) - Era{k + q) 
k,mn 

X E f^"T(fc + q)U:^{k)UJ}l{k)U"^\{k + q), 



al3 



where /„(fc) = 1/(1 + ef'(^"('=-i-''-i) denotes the Fermi- 
distribution function with /i the chemical potential. The 
Hamiltonian (|14l) is diagonalized by a Bogoliubov trans- 
formation yielding the spin-wave dispersion 



u(q) = ^A^q)-B^{q) 



(18) 



of the system in the presence of both, super- and double- 
exchange. 



IV. RESULTS 
A. Classical phase diagram 

As a pre-requisite for the spin-wave expansion, we first 
have to identify the regime where (tt, 0) stripe antifer- 
romagnetism is classically stable. When Jh = or at 




FIG. 2: (Color online) Classical phase diagrams of the degen- 
erate double-exchange model as a function of Jn/ii and J2/J1 
for different filling levels n = 0.05, n — 0.10, and n — 0.15. 
The tight-binding hopping parameters are t2 — — O.lti, = 
0.2ti, and t4 = 0.05ti, the nearest- neighbor superexchange 
Ji = 0.04ti. 



zero filhng (n = 0), (tt, 0) order is stable for J2 > Ji/2. 
In this regime, a finite Hund coupling further stabilizes 
(tt, 0) order since the itinerant electrons are more likely 
to occupy the dyz orbitals which have a larger overlap 
along the ferromagnetic y-direction. Hence the double- 
exchange effectively weakens the spin coupling along the 
y direction and reduces the magnetic frustration. On the 
other hand, for J2 < Ji/2 the (tt, tt) Neel antiferromagnet 
is the classical ground state for Jh — 0. However, from 
the DDEX models for the manganites, it is known^"— 
that even for J2 = 0, a sufficiently strong Hund coupling 
will eventually stabilize (tt, 0) order. 

For J2 < Ji/2, the phase boundary between the two 
classical ground states is determined by the condition 
that the electronic kinetic-energy gain for the orbitally- 
polarized and stripe-ordered configuration equals the 
difference of the magnetic energies, -E(7r.o) ~ ^(tt.tt) = 
2Ji — 4J2- The electronic energies for a given electron 
filling n are easily calculated by diagonalizing the free- 
electron Hamiltonian He pH)) for the two different spin 
configurations. From now on, all the energies will be 
expressed in the unit of ti, the largest of the hopping 
amplitudes. The tight-binding parameters will be chosen 
to be t2 — — O.lii, = 0.2ti, and t4 = 0.05ii through- 
out the paper. We further set Ji = O.OAti so that the 
exchange constants are of the order of 10 meV for a band- 
width of 1 eV, in agreement with both numericaP^- and 
experimental^ observations. 

The resulting phase diagrams for electron fillings n = 
0.05, n — 0.1, and n — 0.15 are shown in Fig. [51 As pre- 
dicted, a sufficiently large Jh can stabihze (tt, 0) stripe 
order even when Neel order is favored by J2 < Ji/2. 
Moreover, the (tt, 0) phase is enhanced for larger filling n 
because the kinetic energy gain increases with the num- 
ber of itinerant electrons. Although similar results have 
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been obtained in the context of the manganites^"— we 
point out that there are crucial differences. Whereas in 
the manganites J2 is neghgible and (tt, 0) antiferromag- 
netism obtains because of the enormously large Hund 
coupling, in the pnictides, Jh is much smaller and of 
the order of the electronic bandwidthj^"— Therefore, the 
large J2, which has been predicted early on^l based on 
the geometry of the ion-arsenic layers, is essential for 
stripe antiferromagnetism in the pnictides. 

In the following calculations we will mostly focus on 
the regime J2 > Ji/2, where (7r,0) order remains clas- 
sically stable for Jh — >• 0. In the calculation, our pri- 
mary interest is in how the spin-wave spectrum is renor- 
malized as we gradually turn on Jh- In principle, when 
J2 < Ji/2, the linear spin- wave calculations can still be 
carried out for a Jh that is large enough to classically 
stabilize the stripe order. 



B. Orbital and spin polarization 

We proceed with a more careful inspection of the elec- 
tronic Hamiltonian T-Lc ([T0| which includes Hund's cou- 
pling to the local moments on a classical level given by 



n^^^ Here we focus on the regime where the classical 



magnetic ground state Ki = exp(jQ-ri) = ±1 is the 
stripe antiferromagnet, Q ~ (tt, 0). Certainly, this is al- 
ways the case for J2 > Ji/2. From the diagonalization 
of 'He pOl) with the same set of hopping parameters used 
previously, we obtain the dispersions of the electronic 
bands £'„(fc) shown in Fig. [3lja) and (b) for Jh = O.lii 
and Jh = l.Oii, respectively. Since the diagonal hopping 
amplitude t4 between different orbitals is non-zero, the 
bands are always momentum-dependent superpositions 
of the two orbitals dyz and dxz as indicated by the color 
coding. For small Jh [Fig. [Sja)], the orbital hybridiza- 
tion is found to be strong for some momenta, for ex- 
ample, along the (0, 0)-(7r, tt) direction. However, in the 
relevant regime where Hund's coupling is of the order of 
the electronic bandwidth, e.g. Jh = IMi [Fig. |3l[b)], the 
gaps between the bands are pronounced and the orbital 
polarization of each band remains almost perfect for all 
momenta. Therefore, from now on, we will denote these 
bands as 'J^z' and 'rfyz' for simplicity. 

We further calculate the total orbital polarization 
"o = HuiPv^,v - Pxz,v) and spin polarization iis = 

Y^a^PaA - Pa,])- Here, pav = {c\ayCiav) IS defined as 
the density of electrons in orbital a = xz,yz with spin 
v =t, 4,. Obviously, these densities sum up to the to- 
tal filling of the bands, n = ^^i^ Pav Note that since 
the densities are defined in the sub-lattice rotated basis, 
v =t corresponds to an electron spin aligned with the 
local moment spin. 

In Fig. ^c), the orbital polarization rio is shown as a 
function of the total filling n for different values of Jh- 
For small n the electrons populate only the lowest band 
with almost perfect dyz character and therefore rio — n. 
Indeed, the slope is very close to one indicating that the 
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FIG. 3: (Color online) (a, b) The dispersions of the itinerant- 
electron bands for Jh — O.lti(a) and Jh = l.Oti(b), along the 
path (0, 0)-(7r, 0)-(7r, 7r)-(0, 0) in the Brillouin zone containing 
one Fe atom per unit cell. The tight-binding parameters are 
chosen to be t2 = —O.lti, ts = Q.2t\, and t4, = O.OSti through- 
out the paper. The bands are colored according to the weight 
of the dyz orbital, (c, d) Total orbital polarization no(c) and 
spin polarization Us (d) as a function of the filling n for various 
Hund's couplings Jh- 



admixture of the d^z orbital to the lowest band is negligi- 
ble and that the orbital polarization can be considered to 
be perfect. The emergent ferro-orbital order, illustrated 
in Fig. [11 is a consequence of the ferromagnetic Hund's 
coupling which wants to align the electron spin with the 
local moments and therefore suppresses the electron mo- 
tion along the antiferromagnetic spin direction. Since 
ti > 1*2 1 the electrons populate the dyz orbitals which 
have a larger overlap along the ferromagnetic y direc- 
tion. Increasing the filling n, the Fermi energy will even- 
tually reach the bottom of the next band with mainly d^z 
character. Above this particular filling n, the orbital po- 
larization is no longer perfect and even starts to decrease 
with n for larger values of Jh- In general, a larger Jh 
increases the energy difference between the two bands, 
and thus increases the filling n up to which the orbital 
polarization is perfect. 

Since Hund's coupling tends to align the spins of the 
itinerant electrons with the local moments, also the spin 
polarization Ug increases with the electron filling n and is 
bigger for larger Jh as shown in Fig. [3Jd). However, for 
small n, the slopes are slightly smaller than one, signaling 
an incomplete spin polarization. This is a consequence 
of the small but finite spin off-diagonal elements in He 
PH)) which do not only depend on the smallest hopping 
amplitude t^, as the orbital hybridization terms, but also 
on t2 and t^. For n > h, the slope of the spin-polarization 
curves is reduced indicating that the spin polarization of 
the dyz orbitals is larger than that of the d^z orbitals. 
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FIG. 4: (Color online) Spin-wave dispersion Ljj{q) in the pres- 
ence of both super- and double-exchange along the path (0, 0)- 
(tt, 0)-(7r, 7r)-(0, 0) for (a) different filling levels n, Jh = 2.0ti, 
and (b) different Hund's couplings Jh, n = 0.1. In both cases, 
Ji = 0.04ti, Ja = 0.6 Ji, and S = 1/2. 



C. Spin- wave spectrum 

Now we set out to calculate the spin-wave spectrum 
by considering the corrections from the itinerant bands 
to the J1-J2 Heisenberg model for the local moments. 
The hopping amplitudes are set to the values used be- 
fore. The spin length S of the local moment is 1/2, 
which reflects the relatively small moment measured by 
the experimenlji^ and is consistent with a local multi- 
plet structure with an orbital degeneracy4S Finally, the 
Heisenberg model is strongly frustrated with Ji = 0.04ii 
and J2 = O.6J1. 

We first calculate the spin-wave spectra for different 
filling levels n with Jh = 2.0ii, shown by Fig. IDJa). In 
this case, the complete orbital polarization is found up to 
n = 0.16 [see Fig.|21[c)]. When n — 0, corresponding to an 
empty itinerant band, our model reduces to an isotropic 
J1-J2 Heisenberg model, where the linear spin- wave ener- 
gies are zero at both, (tt, 0) and (tt, tt). At finite electron 
densities n < n, we observe that the spin-wave energy 
at (7r,7r) is pushed to higher values as n increases. This 
indicates a stabilization of the stripe antiferromagnetism 
over the competing Neel order. We also note a signif- 
icant mode softening along the (0, 0)-(7r, 0) direction at 
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FIG. 5: (Color online) (a) Best fits of the spin-wave spectrum 
of the DDEX for n = 0.1 and Jh ~ S.Oti obtained from a 
spin-only model with effective couplings illustrated as inset, 
(b, c) Effective exchange constants as functions of Hund's 
coupling Jh for bare exchange constants Ji = 0.04ti(b), Ji = 
0.08fi(c), and J2/J1 = 0.6. 



low fillings due to the other finite hopping amplitudes t2 , 
ts, and t^. 

For n > h, not shown in the graph, the spin-wave 
energy at (tt, tt) decreases with n. This behavior re- 
sults from the reduction of the orbital polarization [see 
Fig. EJc)] and hence of the anisotropy, once the itinerant 
electrons populate the next band with mainly dxz char- 
acter. As the filling level continues to rise, the spin-wave 
spectrum becomes unstable suggesting that the system 
may evolve to a different ground state. 

Fig. Hfb) shows the spin- wave spectra for different val- 
ues of the Hund coupling Jh and a filling level oin = 0.1. 
According to Fig. ^c), we have complete orbital po- 
larization for all the Jh's used in Fig. IDJb), except for 
Jh — IMi. As expected, a larger Jh produces stronger 
corrections to the spin- wave dispersion, especially around 
(tt, tt) where the spin- wave energy almost reaches a max- 
imum. In contrast, the dispersion along (0,0) to (7r,0) is 
almost unaffected after Jh reaches a certain value on the 
order of the electronic bandwidth. 



D. Magnetic anisotropy 



In Section llV Cl we saw that the double exchange leads 
to a dramatic change of the magnetic excitation spectra. 
In particular, the spin-wave modes at the Neel antiferro- 
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magnetic wave vector (tt, tt) are almost pushed to a maxi- 
mum consistent with the neutron-scattering data.^^ Since 
orbital ordering leads to a dramatic anisotropy in the 
electronic structure, the ferromagnetic double-exchange 
contribution is expected to be much larger along the y- 
direction along where the ferromagnetic spin chains are 
formed. To quantify the induced magnetic anisotropy, in 
this Section, we fit the spin-wave dispersions calculated 
in the presence of both, super- and double-exchange to 
an effective, anisotropic Heisenberg model. 

To be more specific, in the effective spin-only model, 
nearest-neighbors are coupled by exchanges Jix and Jiy 
along the x and y directions, respectively, and next- 
nearest neighbors by J2 [see inset of Fig. [5l[a)]. Please 
note that different symbols are used here to distinguish 
from the couplings in the original Heisenberg model "Hioc 
([2]). Furthermore, we introduce an additional ferro- 
magnetic exchange J^y < between the third-nearest 
neighbors along the ferromagnetic chains. The utility 
of retaining longer-ranged ferromagnetic couplings has 
been demonstrated in the ferromagnetic double-exchange 
model, where simple nearest-neighbor exchange is un- 
able to reproduce the calculated dispersions from either 
canonical transformations^^'™ or diagrammatic pertur- 
bation theory/^ Also for manganites with so called CE- 
type charge-spin-orbital order, longer-ranged ferromag- 
netic couplings along the ferromagnetic zig-zag chains are 
crucial in order to obtain a good fit of the magnetic ex- 
citation spectral^ 

In the relevant regime Jix > Jiy, Jiy < 2J2, J^y < 
0, where (tt, 0) stripe antiferromagnetism is stable, the 
linear-spin wave dispersion (2; (g) is determined by uj'^iq) = 
A^{q)-B'^{q) with = Ji^ - Jij,(l - cos g^,) 2 J2 - 
"^31/(1 ~ cos 2(7y) and B{q) = Ji^ cos qx + '2J2 cos qx cos qy. 

In Fig. [Sfa), the fitted spin- wave dispersions to the 
spectrum calculated from the DDEX model with n — 0.1 
and Jh = S.Oti are shown. Indeed, the inclusion of 
the longer-ranged ferromagnetic coupling J^y leads to 
a significant improvement of the fit. Moreover, setting 
Jsy = gives an unrealistically large correction to the 
nearest-neighbor exchange along the x-direction {Jix = 
— O.27J1), which has to be compensated by a fairly sub- 
stantial increase in the diagonal exchange, J2 — 0.90Ji. 
This result is certainly unphysical since from the classi- 
cal double-exchange argument, we expect Jix — Ji and 
J2 < J2- In contrast, by including J^y, we obtain more 
physical fitting results. Consequently, in the following, 
all of the fittings will be performed with J^y 7^ 0. 

The resulting effective exchange constants as a func- 
tion of Hund's coupling Jh are shown in Fig. [SJb) and 
(c) for Ji = 0.04ii and Ji = O.OSti, respectively. Note 
that for Ji — 0.04ti, the spectra of the DDEX model are 
shown in Fig. lUb) . In both cases we use n = 0.1, J2/ Ji = 
0.6, and the same tight-binding parameters as before. 
Obviously, for Jh = the effective exchange constants 
have to be identical with the bare ones, Jix = Jij, = Ji, 
J2 = J2, and J3 = 0. Increasing Hund's coupling does 



0.2 




FIG. 6: (Color online) Spin-wave dispersion io{q) in the pres- 
ence of both, super- and double-exchange along the path 
(0, 0)-(7r, 0)-(7r, 7r)-(0, 0) for different Hund's couplings Jh. 
Parameters are the same as in Fig. |4jb), except for J2/J1 = 
0.4. In this regime, (7r,0) order is classically stable for 
Jh > l.Gti. (Inset: the fitted, effective exchange constants 
Jix, Jiy, J2, and Jsy as functions of Jh-) 



not change Jia; , whereas the other exchanges decrease due 
to the different ferromagnetic double-exchange contribu- 
tions. Since ti is the largest hopping amplitude and most 
electrons populate the dyx orbitals, the fastest decrease 
is observed for the coupling Jiy along the ferromagnetic 
chains. We also note that the non-monotonic behavior 
for small Jh is not unphysical, but due to the fact that 
for n — 0.1, complete orbital polarization is not achieved 
until Jh 1.5ii. 

Interestingly, in the case Ji — 0.04ti [see Fig. EJb)], 
the coupling Jiy along the ferromagnetic direction be- 
comes negative for Jh ^ 2.0<i, which completely removes 
the frustration in the effective spin-only model. Such an 
effective negative exchange coupling along the ferromag- 
netic spin direction has been used phenomenologically 
to rationalize the spectra measured by inelastic neutron 
scatteringi^^ Remarkably, around Jh ~ 2.0ti, the ra- 
tios of the three exchange constants Jix, Jiy, and J2, 
agree extremely well with the experimental estimates4^ 
As shown in Fig. [5{c), for Ji — O.OSii, the relative cor- 
rections to the exchange couplings are too small to make 
Jiy ferromagnetic for realistic values of Jh. In order to 
achieve the experimentally observed negative Jiy within 
reasonable parameter space, we require Ji < 0.05ti. In 
fact, inelastic neutron scattering^ suggests that the ex- 
change constants are of the order of 10 meV, which in 
our theory leads to an electronic bandwidth and Hund's 
coupling Jh both of the order of 1 eV, in agreement with 
other experimental estimatesi^ 

Though the parameter regime J2/J1 < 1/2 is most 
likely not relevant to the pnictides, we consider it 
nonetheless for completeness. In this regime, the Hund 
coupling has to exceed a critical value in order to stabi- 
lize (tt, 0) stripe antiferromagnetism. As an example, we 
assume J2 = 0.4Ji and n = 0.1, in which case (tt, 0) order 
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is classically stable for Jh > 1.6ti [see Fig. [2]. The re- 
sulting spin- wave spectra are shown in Fig. [S] for different 
values of Jh- The effective exchange couplings obtained 
from fitting to the spin-only model are shown in the inset. 
Although the behavior is qualitatively similar to the case 
•h/ Ji ~ 0.6, the agreement with the experimental data 
is not as good. Moreover, in order to push the dispersion 
at (tt, tt) close to a maximum, we need a much bigger Jh 
which can not be justified for the pnictides. 

V. DISCUSSION AND CONCLUSION 

In summary, we have studied a DDEX model for the 
iron pnictides which explains the dramatic magnetio^i^ 
and electroni c^*^'^^ anisotropics in these materials. The 
model consists of local moments which are coupled by 
antiferromagnetic superexchanges Ji and J2 between 
nearest- and next-nearest-neighbor spins, respectively, 
and of itinerant electrons in the bands of the degener- 
ate dxz and dyz orbitals. The electrons are coupled to 
the local moments by a ferromagnetic Hund exchange, 
Jh- The system spontaneously develops ferro-orbital or- 
der because of the kinetic energy gained by directing the 
itinerant electrons along ferromagnetic spin chains which 
are stabilized by the double-exchange mechanism. 

Although similar results have been obtained previously 
in the context of the manganites^^"— we point out that 
there is a crucial difference between the two classes of 
materials. Whereas in the pnictides, the local moments 
are inherently quantum {S — 1/2) and Jh is at most of 
the order of the electronic bandwidth, the DDEX models 
for the manganites can be greatly simplified by assum- 
ing an infinitely large Jh and by treating the core spins 
as classical. Because of these approximations, the spin- 
wave excitations of the ferro-orbitally ordered stripe an- 
tiferromagnet have not been addressed in the manganite 
literature. 

In this work, we have explicitly taken the quantum 
nature of the local moments into account and calculated 
the magnetic excitation spectra in the presence of both, 
super- and double-exchange. Over a realistic parame- 
ter range, the calculated spin- wave dispersions are found 
to be in good agreement with the neutron-scattering 
data.— 1^ In particular, we find that the dispersion is 
pushed almost to a local maximum at the competing 
Neel ordering wave vector as seen in the experiment 4^ 
By fitting to an effective spin-only model, we find that 
the coupling along the ferromagnetic direction becomes 
negative, which demonstrates that in the pnictides, the 
double-exchange along the ferromagnetic spin-direction 
can overcompensate the superexchange between the lo- 
cal moments. In this regime, (tt, 0) antiferromagnetism 
is unfrustrated. 



It is feasible that the parent, undoped materials self- 
tune the size of the local moments and the carrier density 
to the point where the (tt, 0) antiferromagnetism is most 
stable. In our theory, this is the case for the optimal fill- 
ing level n = n where the orbital polarization reaches a 
maximum. In fact, the starting Heisenberg model Hioc 
([2|) is likely to be in the regime of strong frustration, 
0.4 < J2/J1 < 0.6, where quantum fluctuations destroy 
long range magnetic order. Only through the interac- 
tion with the itinerant electrons does stripe antiferro- 
magnetism emerge. Electron- or hole-doping the system 
'At n — n diminishes the orbital order and thus increases 
the magnetic frustration, lowering the transition temper- 
ature to a spin-ordered state. 

We conclude by stressing that the degenerate double- 
exchange model motivated and studied in this work qual- 
itatively explains many properties of the undoped and 
slightly doped iron pnictides. The emergent ferro-orbital 
order breaks the in-plane lattice symmetry, thereby driv- 
ing the orthorhombic lattice distortion. Further, or- 
bital order induces a strong electronic anisotropy which 
explains why the structural transition is accompanied 
by dramatic Fermi-surface reconstructions^ and trans- 
port anomalies.— Moreover, recent STM experiments^ 
and in-plane resistivity measurements^!^ have unam- 
biguously demonstrated the spatial anisotropy of the 
electronic state at temperatures below the structural 
transition temperature. Finally, orbital ordering pro- 
duces strong magnetic anisotropy, essential for explain- 
ing the experimentally observed magnetic excitation 
spectra4^ 

During the review of this work, several mean-field 
calculations22r— based on five-orbital Hubbard models 
have shown that the total occupations of the dxz and 
dyz orbitals are close to each other, but that there is a 
significant difference in the density of states for the two 
orbitals at the chemical potential. The ferro-orbital order 
in our theory does not contradict these results. In fact, 
by viewing the itinerant electrons as residual quasipar- 
ticles on top of local moments formed because of strong 
correlations, our model is an effective low-energy theory 
that successfully describes the strongly anisotropic elec- 
tron state around the Fermi surface of the iron pnictides. 
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